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PACS 87. 16. Qp - Pseudopods, lamellipods, cilia, and flagella 
PACS 87. 17. Aa - Modeling, computer simulation of cell processes 
PACS 87. 17. Rt - Cell adhesion and cell mechanics 

Abstract - A cellular bleb grows when a portion of the cell membrane detaches from the under- 
lying cortex under the influence of a cytoplasmic pressure. We develop a quantitative model for 
the growth and dynamics of these objects in a simple two-dimensional setting. In particular, we 
first find the minimum cytoplasmic pressure and minimum unsupported membrane length for a 
stationary bleb to nucleate and grow as a function of the membrane-cortex adhesion. We next 
show how a bleb may travel around the periphery of the cell when the cytoplasmic pressure varies 
in space and time in a prescribed way and find that the traveling speed is governed by the speed 
of the pressure change induced by local cortical contraction while the shape of the traveling bleb 
is governed by the speed of cortical healing. Finally, we relax the assumption that the pressure 
change is prescribed and couple it hydrodynamically to the cortical contraction and membrane 
deformation. By quantifying the phase space of bleb formation and dynamics, our framework 
serves to delineate the range and scope of bleb-associated cell motility and synthesizes a variety 
of experimental observations. 



Introduction — Blebs are protrusions of a cell mem- 
brane driven by local variations in intracellular pressure 
induced by contractility and arc commonly seen in many 
types of cells. Blobbing is a biological analogue of an el- 
ementary mechanical process — the separation of a thin 
membrane from a bulk material. It is thus similar in many 
respects to the formation of a blister in a thin him ad- 
herent to a substrate, but also rather different in that it 
involves a number of active processes associated with ac- 
tive (and regulated) contractile stresses that drive the pro- 
cess, as well as active mechanisms associated with healing 
of these objects via forces at the boundary. Blebs have 
mostly been observed and studied during the processes 
of apoptosis and cytokinesis, but recent studies have re- 
vealed cellular blebbing to be an important mechanism 
contributing to cell motility [l}{5]. For example, migration 
of chemotactic germ cells in some systems is not powered 
by actin polymerization but by cytoplasmic streaming in- 
stead PQ , while some tumor cells form bleb-like protrusions 
to migrate through three-dimensional matrices [5]. It has 
been further suggested that bleb protrusions can cooperate 
with lamellipodium-based protrusions to power cell motil- 



ity [2j[5], and that some types of cells can switch between 
the two modes depending on their environment - 
Thus the focus on actin polymerization-driven protrusions 
in lamellipodia and filopodia on cell motility must be com- 
plemented by considerations of contractility/pressure in- 
duced blebbing before one can decipher the relative con- 
tributions of these modes for whole cell motility. Under- 
standing blebbing quantitatively is a first step in this pro- 
cess. 

Here, we provide a theoretical model for the nucleation, 
expansion, retraction, and large scale movements of blebs. 
A bleb is nucleated when the local hydrostatic pressure 
generated by cortical contraction detaches the cell mem- 
brane from the cortex. Cytoplasmic pressure and flow 
then drives the detached membrane to expand, forming an 
expanding bleb that is not supported by the actin cortex 
except along its edge. This expansion is accommodated 
by further delamination of the cell membrane from the 
cortex, flow of lipid into the bleb through the bleb neck, 
and unwrinkling of excess folded membrane. When the 
bleb expansion slows down, the actin cortex starts to re- 
form underneath the bleb membrane. In non-motile cells, 
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myosin-driven contraction of the reformed cortex retracts 
the bleb, while in motile cells, contraction of the rear of a 
cell leads to a net movement of the cell body towards the 
leading edge where the bleb is formed. In addition, some- 
times an asymmetric reformation of the actin cortex can 
lead to the circus movement of a bleb traveling around the 
periphery of a cell [TU]. During this last process, the actin 
cortex is reformed at one side of the bleb and its contrac- 
tion pushes the cytosol to the other side of the bleb, which 
sometimes leads to further dclamination of the membrane 
at the unconsolidated edge. 

Model — To quantify these phenomena, we focus on a 
simple 2 dimensional model of the dynamics of a cell, con- 
sisting of an active fluid — a contractile cortex bathed in 
a Newtonian fluid, and surrounded by a membrane that 
is detachable from the actin cortex. We assume that the 
membrane is one-dimensional and neglect its natural cur- 
vature, which is reasonable for blebs that are small com- 
pared with the size of the cell. In addition, we focus on 
a single bleb and thus neglect bleb-bleb interactions, as 
well as the adhesion of the bleb to any external scaffolds 
or substrates. For a single bleb, denoting the distance be- 
tween the detached membrane and the cortex as y(x) at 
the position x, the energy per unit width E(x,y,y x ,y xx ) 
of a piece of blistered membrane y{x) is given by 



E 



-^Vxx + -jVx + -jV 



py 



E n dx. 



(1) 



where (.) x = d(.)/dx, and the limits of integration s;. r de- 
fine the left and right boundaries of the bleb, which are 
possibly dynamic. The various terms in the integrand de- 
note, respectively, the bending energy of the membrane 
with flexural rigidity B, the work done by the in-plane 
membrane tension T, the binding energy of the adhesive 
bonds between the membrane and the cortex with spring 
constant k, the work done by the cytoplasmic pressure p, 
and finally the adhesion energy between the membrane 
and the cortex. For consistency, the spring constant k 
must be related to the adhesion energy E a by k = 2E a /l^, 
with l c the maximum length of the membrane-cortex ad- 
hesion bond, with the condition that when y(x) > l c , 
the adhesion bond at x breaks. The corresponding Euler- 
Lagrange equation for the shape of the bleb is then given 
by 



By xxxx - Ty xx + nyH[l - y/Q -p = Q, 



(2) 



where H[a] = l,a < 0;H[a] = 0, a < is related to the 
Heaviside function. In general, the cytoplasmic pressure p 
is non-uniform owing to local variations in the cortical con- 
tractility, and due to the healing of the actin cortex after 
the bleb is formed. We assume that local cortical contrac- 
tion is the main contribution to the change in pressure and 
for now prescribe the pressure of the form of a localized 
traveling pulse given by 



p{x, t) — Ilexp 



(x - v v tf 



(3) 



where II, w^, and x n denote the magnitude, traveling 
speed, and width of the pulse, respectively. This form de- 
couples the relation between actin reformation in the bleb 
to the changes of cortical contractility, a constraint that 
we will relax later on. If v n > 0, the pressure variations 
as well as the bleb boundaries s; jr are time dependent. 
Completing the formulation of the problem for the for- 
mation and dynamics of a bleb, the boundary conditions 
associated with Eq. (2) are 
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for the left boundary and 
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(5) 



for the right boundary where vh- and vh+ are the active 
cortical healing speeds at the left and right boundaries of 
the bleb, controlled by the actomyosin interactions at ei- 
ther boundary. These conditions effectively assume that 
the dominant dissipation mechanisms are associated with 
the the movement of the contact lines, with \i the associ- 
ated dynamic viscosity, and follow by balancing the varia- 
tional derivative of the energy in Eq. (1) with the relative 
velocity of the contact line due to passive mechanical forc- 
ing. Finally, we require y\ S[ r = y x \ Slr =0. 

Our model equations can be made dimensionless by in- 
troducing a characteristic length scale y/ B/T, velocity 
scale T//J,, and energy scale B so that 
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with the dimensionless parameters J, K,y c ,P,x p ,v p ,Vh± 
listed in Table Q] As stated earlier, in general, P = 
P{v P ,Vh±), although we will start by assuming that P 
is not dependent on the contractility or rate of heal- 
ing. Physiologically, the material parameters J, K, and 
y c are properties of adhesion molecules such as the Ezrin, 
Radixin and Moesin (ERM), proteins which are known 
to be key regulators of membrane-cortex interactions and 
signaling jTT]. Thus, modifying the density of these ERM 
proteins in a cell will change the value of J, whereas mod- 
ifying the sequence of the proteins may change their adhe- 
sion strength and rest length and hence K and y c , respec- 
tively. On the other hand, the parameters related to the 
localized traveling pressure pulse, P,x p , and v p are gov- 
erned by the actomyosin contractile activity in the cortex; 
for example, modifying the motor density affects both P 
and x p whereas modifying the actin and actin crosslinkcr 
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density changes v p . Finally, the healing speed Vh± is gov- 
erned by the rate of actin polymerization in the cortex. 
For our calculations, we assume the parameter values as 
listed in Table [IJ chosen to correspond to physiological val- 
ues of B = 10" 19 J [TU], T = 10~ 6 N/m |T0],^ = l(T 2 Pa s, 
E a = lCT 6 N/m, k = 8x 10 10 N/m 3 , l c = 5 x lCT 9 m, II 
ranging from 380 Pa to 475 Pa, x p ranging from 3.2 x 10 -7 
m to 3.2 x 10 -6 m, and the speeds v p and Vh± ranging from 
to 2 x 10" 3 m/s. 

To understand the behavior of the solutions to Eqs. ([6])- 
([5]), we solve them numerically with initial conditions 
s r — si = 2 (or 6.3 x 10~ 7 m, a small patch size) and 
y{x, t = 0) = for all x, using a second-order centered 
finite difference scheme in space and a forward Euler finite 
difference scheme in time. The time-step and the grid-size 
used are dt = 1 x 10~ 4 and dx = 0.025, respectively. 

Results and Discussion — Stationary blebs: We first 
consider the static case when v p = Vh± = 0, and study the 
nuclcation and growth of a bleb as a function of the scaled 
pressure magnitude P and the its spatial width x p . A bleb 
forms when membrane- cortex adhesive bonds are broken 
over a finite length, i.e., y > y c for some x <E (s;,s r ). 
When the pair of values of x p and P are below a thresh- 
old, a bleb cannot form, but if a bleb can be successfully 
nucleated, it will grow into a symmetric stationary bleb. 
The phase diagram of stationary bleb nucleation in P — x p 
space is shown in Fig. [TJ with a phase boundary separat- 
ing values of these parameters that allow for blebs from 
those that do not given by the curve x p ~ (P — P )~ Q , 
in terms of a critical pressure Po and exponent a, which 
depends on the membrane flexural rigidity B. We find 
that a — > 1 when x p S> 1, or equivalently, B/T <C x\ in 
dimensional terms, i.e., when the membrane flexural rigid- 
ity B is small, consistent with prior scaling predictions |10j 
where bending energy was neglected in studying bleb nu- 
cleation. In this limit, the scaling theory [TUj also predicts 
that the critical pressure for bleb nucleation Pq = 2J/y c , 
a relation we also validate numerically, as shown in the in- 
set of Fig. [T] Sample bleb shapes shown in Fig. [T] confirm 
that both the width and the height of a stationary bleb 
increase with P and with x p , as expected. 

Traveling blebs: We now consider the case v p > and 
P > Pq when blebs become non-stationary. For conve- 
nience, we will set Vh+ = and henceforth talk only of 
the relative healing speed Vh = Vh-- A sample numer- 
ical calculation of a traveling bleb is shown in Fig. [U[a- 
b), obtained by solving Eqs. ©-(H]) for parameter values 
that allow for bleb nucleation. We see that a bleb first 
grows and moves in the direction of v p ; its shape eventu- 
ally equilibrates and it travels at a constant speed given by 
v p . However, not all traveling blebs can be sustained. In 
some instances, blebs cannot be formed or retract as they 
travel, even if P > Pq. The traveling bleb decays when ci- 
ther (i) Vh is large, so that the trailing edge catches up with 
the leading edge of the bleb, or (ii) v p is large so that the 
leading edge has a tearing rate that lags behind the loca- 
tion of the pressure pulse. To understand this, we consider 




Fig. 1: Phase diagram of stationary bleb formation when 
v P = Vh± = 0. Critical values of the pressure width x p and the 
pressure magnitude P necessary for bleb formation are found 
to obey x p ~ (P - P )~ a with P = 123.2 and a = 0.59. When 
only large values of x v > 3 are considered (where bending en- 
ergy can be neglected), the exponent approaches unity, namely 
a = 0.73, as previously predicted |10| . The inset shows numer- 
ically calculated Po as a function of adhesion energy J with 
y c kept constant. These numerical results agree well with the 
prediction Po = 2J/y c - Several bleb shapes are also shown. 

the boundary conditions Eqs. ©-([H]) in a frame moving 
with speed v p . We see that increasing v p has the same ef- 
fect as increasing the adhesion energy at the leading edge 
while increasing Vh has the same effect as increasing the 
adhesion energy at the trailing edge. Since a bleb can be 
sustained at equilibrium only when P > P , where Po is 
an increasing function of the adhesion energy </, a high 
value of either v p or Vh will shift Po to a larger value, i.e., 
Po ~ 2{J + v Pi h)/yc, by such an extent that the bleb even- 
tually decays. In other words, there is a critical value of 
the healing speed v p ,h ~ P — 2J/y c above which blebs be- 
come extinct. These results are summarized in Fig. (2Jc), 
where we show two types of non-stationary blebs in the 
3-dimcnsional phase space of P, v p , and Vh, those that are 
steadily traveling and those that decay. 

When a bleb travels with a steady shape and speed, 
our numerical simulations show that this speed is equal to 
v p regardless of the value of Vh, deviating from v p only 
in the initial stage of bleb growth. Furthermore for a 
given traveling speed, the shape of the bleb is determined 
hy Vh] a symmetric bleb exists only when v p = w^/2. 
To understand this, we consider Eqs. ©-(jHl) m a mov- 
ing frame with speed v p . Then, Eq. (JfJ]) becomes time- 
independent and the speeds of the boundaries in the mov- 
ing frame are dsi/dt = — (y X x) 2 /2| S! + J ~ Vp + v h and 
ds r /dt = (y X x) 2 /2|« r — (J + v p ). The steady state so- 
lution of Eq. ((SJ) with the above boundary equations is 
a traveling bleb with speed v p and boundary curvatures 
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lable l: 


Dimensionless parameters in the model described by Eqs. ©-((SI and their values. 


Symbol 


Description 


Dimensional description 


Value 


j 


membrane-cortex adhesion energy 


Ea/T 

kB/T 2 


1 


K 


membrane-cortex spring constant 


5000 


Vc 


critical length of membrane-cortex spring 


lc/V B / T 


0.02 


P 


magnitude of the cytoplasmic pressure 


n^/B/T 3 


100-150 


J, p 


width of the pressure pulse 


xJ^jB/T 


1-10 


v p 


speed of the pressure pulse 




0-20 


Vh± 


cortex-membrane healing speed 


(iv H ±/T 


0-20 









\ 
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Fig. 2: (a) Kymograph showing a traveling bleb obtained from 
Eqs. ([OJ — <[Sj with v p — 1.2 and Vh = 0.6 showing that the trav- 
eling speed is constant; (b) Snapshots of instantaneous bleb 
shapes at t = 1 (dotted), t — 2.5 (dashed), t = 8 (dashed- 
dotted), and t = 15 (solid). The parameter values used were 
P = 130.8, x p = 9.0, Dp = 1.2, Vh = 0.6 and the remaining 
parameter values as listed in Table [Tj (c) Phase diagram de- 
marcating steadily traveling and decaying phases in terms of 
the speed of pressure change v v , speed of cortical healing vh, 
and magnitude of pressure pulse P. Traveling regimes (shaded 
regimes) for P = 132.5, 135, 137.5, and 140 are plotted. The 
traveling speed of a bleb is indicated by the colorbar. The 
solid lines mark the boundaries between the two phases on the 
P planes. When P is increased, the boundary between the 
steadily traveling and decaying phases is shifted outwards. 



(y X x) /2| S! = J-Vp+Vh and (y xx ) /2| f 



J + Vp. There- 
fore, it can be seen that the bleb shape is symmetric when 
u p = Vh/2. 

The role of a nonlocal pressure: In Eq. ([5]), we assumed 
that the localized perturbation for the pressure field is in- 




Fig. 3: Schematic of the model incorporating nonlocal pressure. 
The position of the cell membrane r m is solved from Eqs. (I10II - 
(1131) with the forces due to membrane bending ft, membrane 
tension ft, membrane-cortex adhesion f at j as shown. The inner 
circle of radius r c denotes the cortex contracting with f c , and 
the dashed lines denote membrane-cortex adhesion. 



dependent of other variables. In reality, it is coupled to 
the cortical contraction and healing as well as the shape of 
the cell membrane, as suggested by many experiments [2]. 
To account for this coupling, we consider the interaction of 
the intra- and extra-cellular fluid with the cell membrane 
and solve for the pressure field instead of prescribing it 
as an independent parameter. For relative simplicity, we 
assume that the other parameters associated with local 
cortical contraction and actin reformation, namely v v ,x p , 
and Vh , remain fixed. In the context of this minimal hydro- 
dynamic coupling model, we consider a two-dimensional 
cell immersed in an incompressible viscous fluid, with the 
same viscosity as the cytoplasm, and enclosed by an clas- 
tic membrane uniformly adhered to a permeable but rigid 
actin cortex through soft Hookcan springs; a schematic of 
the model is shown in Fig. [3] 

Initially, the cell membrane and the cortex are as- 
sumed to be circular, with the membrane (cell) given 
by I'm = r r „o? and the rigid cortex given by r c = r c of . 
The membrane-cortex adhesive springs are uniformly dis- 
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tributed and uniformly stretched. This configuration gives 
rise to a uniform positive intracellular pressure. Next, we 
assume that there is a slow global cortical contraction of 
the cell prescribed as r c (t) = r c0 [1 — G (1 — cxp(— t/r))] 
where G is the magnitude of global contraction and 1/r 
characterizes the speed of contraction. This global con- 
traction not only generates a high uniform intracellular 
pressure that assists the local contractile force in nucleat- 
ing and expanding a bleb, but also facilitates bleb volume 
growth. A bleb nucleates when the active cortical contrac- 
tion becomes large enough, modeled as a contractile body 
force 



f c (0, i) = Lexp 



-All 



Vet 



(9) 



of magnitude L and size ~ 1/A, that translates with speed 
v c in the angular direction. 

The presence of a localized body force f c results in a 
non-uniform pressure field p(r) and velocity field u(r) in 
the fluid, which can be deduced from the two-dimensional 
scaled Stokes equation and the equation of continuity 



Vp = 
Vu 



V 2 u 

= 0. 



(10) 

(11) 



We parameterize the cell membrane in the deformed con- 
figuration r m (C) = (x(C),y(C)) T with C = [0,L m ] where 
L m is the perimeter of the membrane. The motion of the 
membrane at position r m is determined from the no-slip 
boundary condition imposed on the Stokes equation, i.e., 



u„ 



dr m 

dt 



(12) 



The permeable rigid cortex that supports the cell mem- 
brane is assumed to be in a position such that at any time 
t, the net force experienced by the cortex is zero. Further- 
more, a spring is assumed to break if the spring energy 
exceeds the membrane-cortex adhesion energy J (i.e. if 
the length of the spring is greater than the critical length 
y c ); this detachment of membrane from the cortex results 
in the nucleation of a bleb. The total body force f (r) in 
Ed.pop is the sum of the forces from membrane bend- 
ing, membrane tension, membrane-cortex adhesion, and 
cytoskeletal contraction, i.e. and is given by 



[f fc (r m ) + f t (r m ) + fad (i"m)] 6(r- r m ) dC 



+ f c (r c )5(r-r c )d(, 
Jo 



(13) 



where L c is the perimeter of the cortex and 5(r) is the 
two-dimensional Dirac delta function. In terms of the 
bending energy density Ef, = 7 2 /2 with curvature 7 = 
{x'y" — y'x"){x 2 + y' 2 )~ 3 / 2 , the stretching energy den- 
sity E t = e 2 /2 with strain e = \J x' 2 + y' 2 — 1, and the 
membrane-cortex spring constant K, the bending force, 



tensile force, and adhesive force can be evaluated in a 
straightforward way. We then use a boundary integral 
method with regularized Stokeslets [12] to solve the result- 
ing equations Eqs. ([T0 | - (fl3|) . while the shape of membrane 
is determined using a forward Euler scheme. 

This allows us to study bleb growth and dynamics by 
following the membrane-cortex detachment as a function 
of the cortical pressure, and track the bleb boundaries un- 
til membrane-cortex adhesion reformation, in accordance 
with Eq. ©. We set the membrane-cortex adhesion ref- 
ormation to start at one boundary after a prescribed time 
interval r r from the time when local cortical contraction 
occurs to prevent membrane-cortex adhesion immediately 
after the initial membrane-cortex detachment and allow 
for bleb growth; r r plays a similar role as the initial patch 
size s r — si\t=o in Eqs. (JT])-© in nucleating the bleb. 

Our enhanced non-local model also leads to steadily 
traveling blebs as shown in Fig. Eta-b). Here, a bleb trav- 
els along with the local contraction speed v c , equivalent 
to the speed of the prescribed pressure pulse motion v p in 
our first linear model. In Fig. Hfc) , we show the region in 
phase space spanned by v c and Vh for which steady trav- 
eling occurs. We find two steady traveling regimes. In 
regime I, shown in Fig.0Jc)), the traveling speed of a bleb 
equals the speed of localized contraction v c and we find 
good agreement in phase space with the linear model as 
shown in Fig. HJc). In regime II, shown in Fig. Etc)), the 
traveling speed depends on the healing speed Vh- This oc- 
curs when the localized contractile body force is relatively 
small but its speed v c is relatively large. For a given rate 
Vh at which actin polymerization and membrane healing 
occurs at the belt edge, the transition from regime I to 
regime II takes place when the contraction speed v c is so 
high such that the tearing rate at the leading edge fails 
to catch up with the local contractile force. This causes 
the bleb to lag behind the cortical contraction so that the 
bleb is then dominated by healing. In a finite domain 
with conserved volume, the combination of the nonlocal 
pressure and active contraction allows such a bleb to be 
sustained via circus movements. However, in an infinite 
domain, these blebs will decay and eventually vanish. If 
v c is increased even further, a bleb cannot be sustained 
even in a closed geometry, and there is a transition from 
the traveling regime II to the extinction phase, as for the 
linear model shown in Fig.^c). 

Conclusion — Our minimal description of the many 
complex phenomena associated with the onset, growth, 
dynamics, and extinction of blebs takes the form of a 
simple theory for the deformations of a thin membrane 
partially attached to a substrate, captured in Eqs. (JBJ- 
ijHJ). We find that bleb growth and bleb motion can be 
modeled by detachment of the cell membrane from the 
underlying actin cortex through the competition between 
the membrane bending energy and the membrane-cortex 
adhesion energy, coupled with a localized pressure term. 
When the pressure pulse is time independent, a bleb will 
grow and equilibrate into a stationary bleb. However, a 
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Fig. 4: (a) Kymograph showing a traveling bleb obtained from 
the Stokes model, with v c = 1.2 and Vh = 2.0; (b) Snapshots of 
instantaneous bleb shapes at t = 2 (dotted), t = 12 (dashed), 
t = 22 (dashed-dotted) and t = 32 (solid). The parameter 
values used were J = 1, K = 20, L = 4, A = 400, r r = 
0.25, v c = 1.2 and v h = 2.0. (c) Phase diagram demarcating 
steadily traveling and decaying phases in terms of the speed of 
localized contraction v c , speed of cortical healing v^, and mag- 
nitude of local contractile force L. Traveling regimes (shaded 
regimes) for L = 4, 4.5, and 5 are plotted. A bleb is successfuly 
formed only when L > 3.35. The traveling speed of a bleb is 
indicated by the colorbar. Two traveling regimes are identi- 
fied: (i) regime in which traveling speed equals to v c (regime 
I characterized by low v c ) and (ii) regime in which traveling 
speed is independent of v c but is dependent on Vh (regime II 
characterized by high v c ). In the vh dependent regime, the 
traveling speed of a bleb is low but increases with increasing 
Vh- As v c further increases, this Vh dependent regime gets nar- 
rower and eventually vanishes at a high enough v c . When L 
is increased, the boundary between the steadily traveling and 
decaying phases is shifted outwards. 

change in cytoplasmic pressure induced by local cortical 
contraction can drive the bleb to move, and lead to active 
blebs that travel steadily. When the pressure pulse associ- 
ated with cortical contraction is an independent function, 
we find our theory is linear and our quantitative results 
are consistent with and complement earlier scaling predic- 
tions [10] . Our theory also allows us to characterize the 
formation and motion of blebs in terms of a phase diagram 
that might be useful in experimental settings. 

More realistically, the actively generated fluid pressure 
is coupled to membrane deformation, local cortical con- 
traction as well the cortical healing dynamics. To un- 



derstand this, we modified our model to account for the 
hydrodynamic coupling between pressure variation and 
membrane deformation. Our simulations show that the 
simple linear model is in agreement with the coupled hy- 
drodynamic model when the localized cytoskcletal con- 
tractile body force is large. In this case local contrac- 
tion dominates healing and the relative contribution to the 
pressure from hydrodynamic causes is small, so that the 
uncoupling of pressure from the other variables is reason- 
able. However, when the contractile forces are relatively 
small, we get active blebbing waves that persist because 
of a combination of passive hydrodynamic and active con- 
tractile and boundary effects, consistent qualitatively with 
prior experiments. Our study sets the stage for bleb driven 
motility, since blebs lead to asymmetric forces so that dy- 
namics of blebs can lead to the motion of whole cells by 
combining the effects described here with the dynamics of 
adhesion of the whole cell. 
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